LORETA With Cortical Constraint: Choosing an Adequate Surface Laplacian Operator

Low resolution electromagnetic tomography (LORETA) is a well-known method for the solution of the l2-based minimization problem for EEG/MEG source reconstruction. LORETA with a volume-based source space is widely used and much effort has been invested in the theory and the application of the method in an experimental context. However, it is especially interesting to use anatomical prior knowledge and constrain the LORETA's solution to the cortical surface. This strongly reduces the number of unknowns in the inverse approach. Unlike the Laplace operator in the volume case with a rectangular and regular grid, the mesh is triangulated and highly irregular in the surface case. Thus, it is not trivial to choose or construct a Laplace operator (termed Laplace-Beltrami operator when applied to surfaces) that has the desired properties and takes into account the geometry of the mesh. In this paper, the basic methodology behind cortical LORETA is discussed and the method is applied for source reconstruction of simulated data using different Laplace-Beltrami operators in the smoothing term. The results achieved with the different operators are compared with respect to their accuracy using various measures. Conclusions about the choice of an appropriate operator are deduced from the results

Ongoing EEG artifact correction using blind source separation

Objective: Analysis of the electroencephalogram (EEG) for epileptic spike and seizure detection or brain-computer interfaces can be severely hampered by the presence of artifacts. The aim of this study is to describe and evaluate a fast automatic algorithm for ongoing correction of artifacts in continuous EEG recordings, which can be applied offline and online. Methods: The automatic algorithm for ongoing correction of artifacts is based on fast blind source separation. It uses a sliding window technique with overlapping epochs and features in the spatial, temporal and frequency domain to detect and correct ocular, cardiac, muscle and powerline artifacts. Results: The approach was validated in an independent evaluation study on publicly available continuous EEG data with 2035 marked artifacts. Validation confirmed that 88% of the artifacts could be removed successfully (ocular: 81%, cardiac: 84%, muscle: 98%, powerline: 100%). It outperformed state-of-the-art algorithms both in terms of artifact reduction rates and computation time. Conclusions: Fast ongoing artifact correction successfully removed a good proportion of artifacts, while preserving most of the EEG signals. Significance: The presented algorithm may be useful for ongoing correction of artifacts, e.g., in online systems for epileptic spike and seizure detection or brain-computer interfaces.Comment: 16 pages, 4 figures, 3 table

Sensitivity of a 29-Channel MEG Source Montage

Publisher Copyright: © 2022 by the authors. Licensee MDPI, Basel, Switzerland.In this paper, we study the performance of a source montage corresponding to 29 brain regions reconstructed from whole-head magnetoencephalographic (MEG) recordings, with the aim of facilitating the review of MEG data containing epileptiform discharges. Test data were obtained by superposing simulated signals from 100-nAm dipolar sources to a resting state MEG recording from a healthy subject. Simulated sources were placed systematically to different cortical locations for defining the optimal regularization for the source montage reconstruction and for assessing the detectability of the source activity from the 29-channel MEG source montage. The signal-to-noise ratio (SNR), computed for each source from the sensor-level and source-montage signals, was used as the evaluation parameter. Without regularization, the SNR from the simulated sources was larger in the sensor-level signals than in the source montage reconstructions. Setting the regularization to 2% increased the source montage SNR to the same level as the sensor-level SNR, improving the detectability of the simulated events from the source montage reconstruction. Sources producing a SNR of at least 15 dB were visually detectable from the source-montage signals. Such sources are located closer than about 75 mm from the MEG sensors, in practice covering all areas in the grey matter. The 29-channel source montage creates more focal signals compared to the sensor space and can significantly shorten the detection time of epileptiform MEG discharges for focus localization.Peer reviewe

The Discontinuous Galerkin Finite Element Method for Solving the MEG and the Combined MEG/EEG Forward Problem

In Electro- (EEG) and Magnetoencephalography (MEG), one important requirement of source reconstruction is the forward model. The continuous Galerkin finite element method (CG-FEM) has become one of the dominant approaches for solving the forward problem over the last decades. Recently, a discontinuous Galerkin FEM (DG-FEM) EEG forward approach has been proposed as an alternative to CG-FEM (Engwer et al., 2017). It was shown that DG-FEM preserves the property of conservation of charge and that it can, in certain situations such as the so-called skull leakages, be superior to the standard CG-FEM approach. In this paper, we developed, implemented, and evaluated two DG-FEM approaches for the MEG forward problem, namely a conservative and a non-conservative one. The subtraction approach was used as source model. The validation and evaluation work was done in statistical investigations in multi-layer homogeneous sphere models, where an analytic solution exists, and in a six-compartment realistically shaped head volume conductor model. In agreement with the theory, the conservative DG-FEM approach was found to be superior to the non-conservative DG-FEM implementation. This approach also showed convergence with increasing resolution of the hexahedral meshes. While in the EEG case, in presence of skull leakages, DG-FEM outperformed CG-FEM, in MEG, DG-FEM achieved similar numerical errors as the CG-FEM approach, i.e., skull leakages do not play a role for the MEG modality. In particular, for the finest mesh resolution of 1 mm sources with a distance of 1.59 mm from the brain-CSF surface, DG-FEM yielded mean topographical errors (relative difference measure, RDM%) of 1.5% and mean magnitude errors (MAG%) of 0.1% for the magnetic field. However, if the goal is a combined source analysis of EEG and MEG data, then it is highly desirable to employ the same forward model for both EEG and MEG data. Based on these results, we conclude that the newly presented conservative DG-FEM can at least complement and in some scenarios even outperform the established CG-FEM approaches in EEG or combined MEG/EEG source analysis scenarios, which motivates a further evaluation of DG-FEM for applications in bioelectromagnetism

Relative Source Power: A novel method for localizing epileptiform EEG discharges

Scherg M, Schulz R, Berg P, et al. Relative Source Power: A novel method for localizing epileptiform EEG discharges. Clinical Neurophysiology . 2021;133:9-19.OBJECTIVE: To validate relative source power (RSP) imaging of extratemporal interictal epileptiform discharges (IEDs).; METHODS: The accuracy of RSP was validated in a cohort of patients with extratemporal focal epilepsy and a confined epileptogenic lesion (<19cm3) using distance to the lesion, concordance with resected area and postoperative outcome. Performance was compared with three conventional methods: voltage maps, equivalent current dipole and a distributed source model.; RESULTS: Thirty-three of 41 consecutive patients (80%) had IED averages suitable for analysis. While the peak negativity in voltage maps localized above the epileptogenic lesion only in 18 cases, RSP-maps matched in 29 cases (88%, p<0.0026). Source localization showed a median distance of 9.8mm from the lesion. Source-regions with 20mm radius included 98% of all source-to-lesion distances. In the 21 surgical cases, outcome showed a sensitivity of 82.35% and specificity of 50% without significant differences between the three source imaging methods.; CONCLUSIONS: RSP-maps provide a rapid, intuitive and more accurate source estimation than voltage maps. At sublobar level, RSP localizes with an accuracy similar to conventional methods and results of previous studies.; SIGNIFICANCE: The definition of a source region with 20mm radius helps in guiding further exploration in extratemporal focal epilepsy. Copyright © 2021 International Federation of Clinical Neurophysiology. Published by Elsevier B.V. All rights reserved

The Discontinuous Galerkin Finite Element Method for Solving the MEG and the Combined MEG/EEG Forward Problem

In Electro- (EEG) and Magnetoencephalography (MEG), one important requirement of source reconstruction is the forward model. The continuous Galerkin finite element method (CG-FEM) has become one of the dominant approaches for solving the forward problem over the last decades. Recently, a discontinuous Galerkin FEM (DG-FEM) EEG forward approach has been proposed as an alternative to CG-FEM (Engwer et al., 2017). It was shown that DG-FEM preserves the property of conservation of charge and that it can, in certain situations such as the so-called skull leakages, be superior to the standard CG-FEM approach. In this paper, we developed, implemented, and evaluated two DG-FEM approaches for the MEG forward problem, namely a conservative and a non-conservative one. The subtraction approach was used as source model. The validation and evaluation work was done in statistical investigations in multi-layer homogeneous sphere models, where an analytic solution exists, and in a six-compartment realistically shaped head volume conductor model. In agreement with the theory, the conservative DG-FEM approach was found to be superior to the non-conservative DG-FEM implementation. This approach also showed convergence with increasing resolution of the hexahedral meshes. While in the EEG case, in presence of skull leakages, DG-FEM outperformed CG-FEM, in MEG, DG-FEM achieved similar numerical errors as the CG-FEM approach, i.e., skull leakages do not play a role for the MEG modality. In particular, for the finest mesh resolution of 1 mm sources with a distance of 1.59 mm from the brain-CSF surface, DG-FEM yielded mean topographical errors (relative difference measure, RDM%) of 1.5% and mean magnitude errors (MAG%) of 0.1% for the magnetic field. However, if the goal is a combined source analysis of EEG and MEG data, then it is highly desirable to employ the same forward model for both EEG and MEG data. Based on these results, we conclude that the newly presented conservative DG-FEM can at least complement and in some scenarios even outperform the established CG-FEM approaches in EEG or combined MEG/EEG source analysis scenarios, which motivates a further evaluation of DG-FEM for applications in bioelectromagnetism